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We address the stability of the surface phases that occur on the C-side of 3C-SiC(lll) at the 
onset of graphene formation. In this growth range, experimental reports reveal a coexistence of 
several surface phases. This coexistence can be explained by a Si-rich model for the unknown (3x3) 
reconstruction, the known (2x2)c adatom phase, and the graphene covered (2x2)c phase. By 
constructing an ah initio surface phase diagram using a van der Waals corrected density functional, 
we show that the formation of a well defined interface structure like the “buffer-layer” on the Si side 
is blocked by Si-rich surface reconstructions. 


Graphene grown on silicon carbide (SiC) is one of 
the most promising material combinations for future 
graphene applications.[1-6] On SiC, graphene growth is 
achieved by thermal decomposition of the substrate. [7- 
9] The electronic properties of few-layer graphene films 
grown on the C-side of the polar SiC surface are simi¬ 
lar to those of an isolated monolayer graphene film with 
very high electron mobilities. [10-12] However, control¬ 
ling the layer thickness of the graphene films remains a 
challenge. [13] While some groups report the successful 
growth of large-scale monolayer graphene [14, 15], other 
reports suggest that the pure monolayer growth regime is 
difficult to achieve on the C face of SiC. [13] This is very 
different from the Si side, where nearly perfect, mono- 
layer graphene films can be grown over large areas. [2, 5] 
Investigating the relative phase stability of the competing 
surface phases in the thermodynamic range of graphiti- 
sation is an important step for a better understanding of 
graphene growth. 

For epitaxial graphene films on the Si-side of 3C- 
SiC(lll), we have recently shown by ah initio atom¬ 
istic thermodynamics that individual phases, the (6v^ x 
6V3)-R30° zero layer and monolayer graphene (ZLG and 
MLG) can form as near equilibrium phases under certain 
external conditions (represented by the C and Si chem¬ 
ical potentials that are controlled by temperature and 
background gas pressure in experiment, see also Eq. (1) 
below).[16] What is not a priori clear is whether on the 
C face of SiC, graphene films can also be thermodynam¬ 
ically stable. To address this question, atomistic models 
are required for the different competing phases, in partic¬ 
ular for the C rich conditions close to the graphitization 
regime. 

We here present first-principles evidence that the for¬ 
mation of monolayer graphene films on the C face is hin¬ 
dered by stable Si rich phases. This is a major and un¬ 
expected difference to the case of the Si face, where the 
formation of a C-rich, so-called “buffer layer” phase is 
actually aided by the formation of heterogeneous C-Si 


bonds. To shed light on the phase mixture at the graphi- 
tisation limit on the C face, we use a possible model of 
the unknown (3x3) reconstruction. The central feature 
of this model is a capping layer of Si atoms, minimiz¬ 
ing dangling bonds in the same way as the known Si-rich 
(3x3) phase on the Si face of SiC.[17, 18] 

On the C-side, a series of different surface structures 
have been observed during annealing. [7, 19-27] Graphene 
growth starts either with a Si rich (2x2) phase in an Si 
rich environment[21, 22] or with an oxidic (a/3 x ^/3) 
reconstruction. [21, 28] In the absence of a Si back¬ 
ground like disilane, a (1x1) phase is observed, which 
exhibits the periodicity of bulk SiC underneath a disor¬ 
dered oxidic layer. [19] Continued heating leads to a (3x3) 
phase. Using qualitative low-energy electron diffraction 
(FEED), Bernhardt et al. [22] showed that the (3x3) 
reconstructions originating from different starting struc¬ 
tures and environments are equivalent. [22] Eurther an¬ 
nealing leads to a (2x2) Si ad-atom phase, referred to as 
(2x2)c (notation taken from [22]). Just before graphene 
forms on the surface, a coexistence of the two surface 
phases, the (3x3) as well as the (2x2)c, is observed. 
[22, 29] In addition, different groups reported strong ex¬ 
perimental evidence that both reconstructions persist un¬ 
derneath the graphene films with the (3x3) phase gradu¬ 
ally fading, but never disappearing. [26, 27, 29-34] While 
the atomic structure of the (2x2)c reconstruction was 
resolved by quantitative FEED [23], the (3x3) recon¬ 
struction and the graphene/ SiC interface remain a puz¬ 
zle. In the past, several structural models were suggested 
for the (3x3) phase on the C face. [20, 21, 27, 35] We 
here summarise its experimentally observed characteris¬ 
tics. Its stoichiometry was found to be Si-rich by Auger 
electron spectroscopy (AES).[20, 22] The corresponding 
filled state scanning tunneling microscope (STM) image 
is consistent with three adatoms residing at the same 
height. [27, 35] Scanning tunnelling spectroscopy shows a 
semiconducting surface with a 1.5 eV band gap. [25] 

Eor the graphene/SiC interfaces on the C-side two dif- 
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FIG. 1. Comparison of the surface energies relative to the 
bulk terminated (1x1) phase as a function of the C chemical 
potential within the allowed ranges (given by diamond Si and 
graphite C). Shaded areas indicate chemical potential values 
outside the strict thermodynamic stability limits. The surface 
energy diagram includes structure models proposed earlier for 
the C face [(b) [23] (d)[21], (e) and (h)[20], (f)[27], (g)[35]] and 
models adapted from the Si face [(a) [17] and (c)[36]]. 


ferent scenarios have been invoked, a) The first carbon 
layer is strongly bound to the substrate. [13, 37, 38] In 
this scenario, the Si sublimation rate during graphene 
growth is controlled by either working in an inert gas at¬ 
mosphere [2], by using a confined geometry [33], or by 
providing an external Si gas phase [ 8 , 38] for example 
disilane. It is not clear, however, if the different groups 
observe the same structures. Based on their FEED 
data, Srivastava et al. [38] proposed a (\/43 x \/43)- 
R±7.6° SiC substrate with a ( 8 x 8 ) carbon mesh rotated 
by 7.6° with respect to the substrate (’\/43-R7.6°’).[38] 
The detailed atomic structure of this interface is not 
known. [39] b) Eor samples prepared under ultra-high vac¬ 
uum conditions, the first carbon layer is weakly bound 
to the substrate, showing the characteristic behavior of 
the TT-band at the iF-point of the Brillouin zone. Here 
an inhomogeneous interface is present since the (3x3) 
as well as a ( 2 x 2 )c reconstruction is observed under¬ 
neath the graphene layers. [26, 31, 32] A recent study on 
graphene grown by molecular beam epitaxy (MBE) on 
the C-side exhibited the same structural characteristics 
as graphene grown by high-temperature annealing. [34] 
This is a strong indication that indeed the (2x2)c as well 
as the (3x3) reconstruction prevails below the graphene 
films. 

To our knowledge, the different (3x3) models sug¬ 
gested in the literature and likewise the remaining phases 
including the graphene/SiC interfaces have not yet been 
placed in the context of a surface phase diagram. We 
employ density-functional theory (DET) using the EHI- 
aims all-electron code [40, 41] with the ELBA eigensolver 
library. [42, 43] We use the van der Waals (vdW) corrected 
[44] Perdew-Burke-Enzerhof (PBE) generalised gradient 


approximation[45] (PBE+vdW) and the Heyd-Scuseria- 
Ernzerhof hybrid functional (HSE06-f-vdW) [44, 46] for 
the exchange correlation functional. Unless otherwise 
noted the calculations are non-spinpolarised. Technical 
parameters and bulk lattice constants are listed in the 
supplemental material (SM)[47]. In contrast to hexago¬ 
nal polytypes, at the surface of 3C-SiC only one type of 
stacking order is present [19]. Although the growth pro¬ 
cess [48] and the electronic structure [49] differ among 
polytypes, the surface reconstruction does not seem to 
be affected [18, 49, 50]. 

A good indicator for finding the most likely (3x3) re¬ 
construction or interface structures is a comparison of 
the respective surface free energies as formulated in the 
ab initio atomistic thermodynamics approach. [51-55] We 
neglect vibrational and configurational entropy contri¬ 
butions to the free energy, although in the coexistence 
region they might lead to small shifts. In the limit of 
sufficiently thick slabs, the surface energy 7 of a two- 
dimensional periodic SiC slab with a C face and a Si face 
is given as 


7Si face + 7C face = ^ - A'si/^Si “ ^CMc) • (1) 

A'si and Nq denote the number of Si and C atoms in 
the slab, respectively, /isi and fic refer to the chemical 
potentials of Si and C. The stability of the SiC bulk dic¬ 
tates /isi+Mc = ^siC^ where is the total energy of a 

bulk SiC unit cell. All surface energies are given in eV per 
area {A) of a ( 1 x 1 ) SiC unit cell. The letter E denotes 
total energies for a given atomic geometry throughout 
this work. The chemical potential limits of the C and 
Si reservoirs are fixed by the requirement that the un¬ 
derlying SiC bulk is stable against decomposition [16], 
leading to: — E^A^ < Me < Because of 

the close competition between the diamond and graphite 
structure for C [56-58], we include both limiting phases 
in our analysis. 

The surface energies of the (2x2)c surface model by 
Seubert et al. and the different models for the SiC- 
(3x3) reconstruction are shown as a function of A/ic = 
Me — Eq^^^ in Eig. 1. The structure with the lowest en¬ 
ergy for a given Amc corresponds to the most stable 
phase. We here briefly discuss the alternative surface 
phases of the SiC C face (for more details see SM [47]). 
Hoster, Kulakov, and Bullemer suggested a geometric 
configuration for the (3x3) reconstruction on the basis 
of STM measurements without specifying the chemical 
composition. [20] Highest in surface energy is the varia¬ 
tion suggested by Hiebel et al. [27] labeled followed 
by a carbon rich composition suggested by Deretzis and 
La Magna [35] labeled g. We added a modification with 
all adatoms chosen to be Si, labeled e. Of all the Hoster- 
type models that we tested, this is the most stable chem¬ 
ical composition. Hiebel et al. suggested a new model, 
labeled /. [27] Li and Tsong proposed a tetrahedrally 
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FIG. 2. I: The a) and b) phases from Fig. 1 calculated using 
the HSE06 exchange-correlation functional with fully relaxed 
structures and unit cells. II: The geometry of the Si twist 
model and simulated constant current STM images for occu¬ 
pied and empty state of the Si twist model (unit cell shown in 
red). The three points of interest (A, B, C) marked by arrows 
are labeled according to Hiebel et a/. [27]. Ill: Spin-polarized 
density of states. 


shaped cluster as reconstruction. [21] They suggested a Si 
and C rich configuration. We tested both configurations 
and included an additional Si tetrahedron (for details see 
SM [47]). Here, we include only the most stable cluster 
formed by 4 Si atoms, labeled d. Finally, we added a 
model originally proposed as a Si rich structure for the 
6H-SiC(0001)-(3x3) reconstruction by Kulakov, Henn, 
and Bullemer [36], labeled c. To summarise Fig. 1, all 
the alternative models we tested are too high in energy 
at the graphite line to coexist with the (2x2)c adatom 
model. 

We next show that a conceptual model for the (3x3) 
reconstruction based on a Si-rich termination performs 
much better at explaining the various surface character¬ 
istics. To create a plausible termination, we base our 
model on the Si twist model [17, 18], known from the 
3C-SiC(lll)-(3x3) reconstruction. In Fig. 2 (II) its ge¬ 
ometry is shown in a side view and from atop. The top 
bulk C layer is covered by a Si adlayer forming hetero¬ 
geneous Si-C bonds. Three Si adatoms form a triangle 


twisted by 7.7° with respect to the top SiC layer. In com¬ 
parison, the twist angle on the Si side amounts to 9.3°. 
The topmost Si adatom is positioned on top of the trian¬ 
gle. In the surface diagram (Fig. 1), this phase has the 
lowest energy of all previously proposed (3x3) models. 
Its formation energy crosses that of the (2x2)c phase 
just at the C rich limit (graphite) of the chemical poten¬ 
tial. To coexist with the (2x2)c phase and to be present 
at the onset of graphite formation, the (3x3) phase has 
to cross the graphite line very close to the crossing point 
between the graphite line and the (2x2)c phase. The Si 
twist model shown in Fig. 2 satisfies this condition. 

In agreement with the AES experiments our model is 
Si-rich. We also compared the surface energetics of the 
coexisting phases (2x2)c and the Si twist model using 
the higher level HSE06+vdW hybrid functional with fully 
relaxed structures and unit cells, shown in Fig. 2 (I). As 
can be seen, the phase coexistence does not depend on 
the chosen functional. However, to distinguish between 
a phase coexistence or a close competition between the 
two phases, the inclusion of entropy terms would be the 
next step. [59] 

In Fig. 2 (HI), we show the spin-polarised electronic 
density of states (DOS) for the spin down (in red) and 
spin up (in blue) channel. While the spin down surface 
state gives rise to a peak in the band gap above the Fermi- 
level, the spin up surface state is in resonance with the 
SiC bulk states. The DOS clearly demonstrates that the 
surface is semiconducting in agreement with experiment, 
featuring a band gap of 1.12 eV. Thus, the conceptual Si 
twist model - inspired by the Si side - appears to satisfy 
the existing experimental constraints well. 

Furthermore, our simulated STM images Fig. 2 (H) 
reproduce the measured height modulation, [27] but the 
experimentally observed difference in intensity between 
occupied and empty state images is not captured by our 
simulated images. The disagreement in the STM images 
might be an indication that a different structure is ob¬ 
served in the STM measurements. We started an exhaus¬ 
tive structure search [60] to find a surface model that is 
even lower in energy than the Si twist model and that 
reproduces all experimental observations, including the 
STM images. However, if an alternative model were to 
be found, its surface energy would have to be close to the 
Si twist model at the graphite line to still coexist with the 
(2x2)c phase. As a result, it would very likely coexist 
with the Si twist model. 

In the following we use the Si twist model as a represen¬ 
tative model to shed light on the SiC-graphene interface 
on the C face. In particular, we will make a simple qual¬ 
itative argument why the C-rich ZLG interface, critical 
towards MLG formation on the Si face, does not form. 
Figure 3 shows the PBE+vdW surface energies of four 
different interface structures, the (2x2)c surface phase 
and the proposed (3x3) Si twist model 

As a first step we constructed an interface structure 
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|c) 6\/3-(2x2)c interf. | |d) 6\/3-(3x3) interf. | 



|e) 6-/3-R30° interf! | |f) -/43-R7.6° inter£ | 




FIG. 3. Comparison of the surface energies for four differ¬ 
ent interface structures of 3C-SiC(lll), relative to the bulk- 
truncated (1x1) phase, as a function of the C chemical poten¬ 
tial within the allowed ranges (given by diamond Si, diamond 
C or graphite C, respectively), using the graphite limit as zero 
reference. In addition, the known (2x2)c reconstruction and 
the (3x3) Si twist model are shown. 

similar to the ZLG phase known from the Si face - a 
6 v/3-R30°-interface, labeled e in Fig. 3. This structure 
crosses the graphite line 0.46 eV above the crossing point 
of the (3x3) Si twist models rendering it unstable. As a 
second structure, we included a purely C based model of 
the v^iS-R 7.6°-interface, labeled / in the surface phase 
diagram of Fig. 3. In our calculations, this structure is 
even higher in energy than the 6\/3-R30°-interface. 

The 6 \/3-R30°- and \/43-R 7.6°-interface are models 
of a strongly bound interface. Since the Gibbs energy of 
formation for SiC is quite large (—0.77 eV in experiment 
[61], —0.56 eV in DFT-PBE-fvdW and —0.59 eV in DFT- 
HSE06+vdW), the formation of SiC bonds is favorable. 
This explains why the (2x2)c and (3x3) Si twist model 
are more stable, because they contain a large number of 
Si-C bonds. Conversely, the 6v/3-R30°- and the y^iS-R 
7.6°-interface are made up of energetically less favorable 
C-C bonds, which increases the surface energy consider¬ 
ably. 

Eor the weakly bound interface a (2x2) and (3x3) 
FEED pattern was observed underneath graphene [26, 
31, 32]. A typical feature of the FEED structure is a 
ring like pattern originating from the rotational disorder 
of graphene films grown on the C face. [10, 29, 31, 62] 
To model the interface, we limited our study to a 30° 
rotation between the substrate and the graphene film. 
This choice was motivated by the FEED study of Hass 
et a/. [10], who showed that graphene sheets on the C 


face appear mainly with a 30° and a ±2.2° rotation. 
A 30° rotation has also be seen in STM measurements 
for the graphene covered (2x2) and (3x3) phases [31], 
from here on called (2x2)g and (3x3)g- We therefore 
chose a (6\/3 x 6\/3) SiC supercell covered by a (13 x 13) 
graphene cell rotated by 30° with respect to the sub¬ 
strate. 

The (2x2)G-mterface covers 27 unit cells of the (2x2)g 
reconstruction (labeled c in Eig. 3). It crosses the (2x2)g 
reconstruction just to the right of the graphite limit at 
a chemical potential of 2 meV and a surface energy of 
—0.69 eV. This finding demonstrates that the observed 
(2x2) FEED pattern underneath the graphene layer is 
indeed consistent with the well known (2x2)g reconstruc¬ 
tion. Our model of the graphene covered (3x3) phase 
consist of the same (13x13) graphene supercell, covering 
12 units of the (3x3) Si twist model (labeled d in Eig. 3). 
The surface energy difference between (2x2)g- and the 
(3x3)G-interface amounts to 0.13 eV at the graphite line, 
favoring the (2x2)G-interface. 

To put our results into context, we revisit the growth 
process. In experiment, growth starts with a clean 
(3x3) reconstruction. The sample is annealed until 
the surface is covered by graphene. At this stage, the 
(3x3) reconstruction is the dominant phase underneath 
graphene. [27] However, a shift from the SiC (3x3) to 
the (2x2)g surface reconstruction at the graphene/ SiC 
interface can be stimulated by an additional anneal¬ 
ing step at a temperature below graphitisation (950C° 
- 1000C°) leaving the graphene layer unaffected.[27] 
Graphene growth starts at a point where bulk SiC de¬ 
composes. The interface is determined by the momen¬ 
tary stoichiometry at which the sublimation stopped - 
a coexistence of different phases is observed. The fi¬ 
nal annealing step shifts the chemical potential into a 
regime in which SiC bulk decomposition stops and the 
graphene layer does not disintegrate, but the SiC sur¬ 
face at the interface moves closer to local equilibrium, 
in agreement with the phase diagram in Eig. 3, forming 
an interface structure between the (2x2)g reconstruction 
and graphene. 

In summary, we shed new light on the central aspects 
of the thermodynamically stable phases that govern the 
onset of graphene formation on SiC (111). In our view, 
the main difference between the Si face and the C face is 
the fact that Si-terminated phases are more stable on the 
C face due to the formation of heterogeneous Si-C bonds. 
The Si rich phases are thus stable practically up to the 
graphite line, allowing the surface to form graphene only 
at the point where bulk SiC itself is no longer stable. 
This makes monolayer graphene growth on the C side 
more difficult than on the Si face. 
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supported by the DEG collaborative research project 951 
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